Topological Superfluid Transition Induced by Periodically Driven Optical Lattice 
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We propose a scenario to create topological superfluid in a periodically driven two-dimensional 
square optical lattice. We study the phase diagram of a spin-orbit coupled s-wave pairing superfluid 
in a periodically driven two-dimensional square optical lattice. We find that a phase transition from 
a trivial superfluid to a topological superfluid occurs when the potentials of the optical lattices are 
periodically changed. The topological phase is called Floquet topological superfluid and can host 
Majorana fermions. 
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I. INTRODUCTION 

Optical lattice system has gradually become a promis- 
ing platform to study many-body quantum systems be- 
cause of lots of significant advances in cold-atom exper- 
iments [1]. In particular, recent theoretical and experi- 
mental progress in laser-induced-gauge-field \2-i7\ makes 
it a hot spot to study topological quantum states in 
cold atoms system [Sl-fllj. Recently, topological quan- 
tum states have attracted considerable interest in con- 
densed physics; however, subject to the compounds' nat- 
ural properties [l2|, we have to rely on serendipity in 
loo king for topological materials in solid-state structures 
|13l - [la |. In contrast, one can engineer the Hamiltonian 
of an optical lattice system to realize variant quantum 
phase states I17H19J . 

In this paper, we show that periodically driven per- 
turbations may give rise to a phase transition from a 
trivial superfluid to a topological one, which carries the 
hallmark with topological protected gapless edges on 
the boundaries of the system. Time-periodic depen- 
dent Hamiltonian can be described by Floquet's theo- 
rem, which is used to explain quantized adiabatic pump- 
ing phenomena [2fl-|2j]. Recently, it demonstrated that 
the phase transition from a superfluid to a Mott insula- 
tor in one-dimensional Bose-Hubbard model can be in- 
duced by a periodically driven optical lattice [2J]. We 
extend this phase transition mechanism to explore the 
topological phase transition in a two-dimensional opti- 
cal lattice. We study the phase diagram of a spin-orbit 
coupled s-wave pairing superfluid in a periodically driven 
two-dimensional (2D) square optical lattice. We find that 
a topological phase transition from a trivial superfluid 
to a topological superfluid can be induced in periodi- 
cally modulated optical lattices. The topological phase 
is called Floquet topological superfluid [25-29] and can 
host Floquet Majorana fermions. It was proposed that a 
topological phase can be realized in a BCS s-wave super- 



fluid of ultracold fcrmionic atoms in the presence of both 
a Rashba spin-orbit ( SO) interaction and a large perpen- 
dicular Zeeman field |30l - l33j : however, the Rashba spin- 
orbit coupling and a large perpendicular Zeeman field 
are hard to be simultaneously realized for cold fermionic 
atoms [32, |33| . We will prove that if one replaces the Zee- 
man field by a periodically driven optical lattice [2^, a 
spin-orbit coupled BCS s-wave superfluid will still allow 
a realization of topological superfluid through modifying 
the oscillating amplitude (or modulation strength) of op- 
tical lattice. Therefore, we provide an alternative method 
to create an important topological superfluid which can 
host Majorana fermions. 

The paper is organized as follows: In Sec. II, we in- 
troduce the s-wave superfluid model in a square optical 
lattice in the presence of both a Rashba SO coupling 
and a periodically modulated optical lattice potential. 
A Zeeman-magnetic-flcld-like term will be derived un- 
der the first order approximation; In Sec. Ill we present 
a two-band approximation and explain the topological 
phase transition at the F point in the first Brillouin zone 
(BZ). At last, we give a brief summary in Sec. IV. 



II. MODEL 

The tight-binding Hamiltonian, which describes an s- 
wave superfluid of neutral fermionic atoms in a 2D optical 
square lattice, is given by 



where 



H (t) = Ho + Hd it) 



Ho = -tY, c\cj - iA ^ c\ (a X d, 



(1) 



(2) 



and 



tions with the scalar product given by 



Hd{t) = K{t)- > r,:cic,: 



E' 



(3) 



((•|-»-^^ dti 



(6) 



Here t is the hopping amphtude between the nearest 
neighbor hnk {i,j), cl^(cl^,clA with c|„ (da) denot- 
ing the creation (annihilation) operator of a fermionic 
atom with pseudospin a (up or down) on lattice site i. 
The second term in Eq. ([5]) represents a Rashba SO cou- 
pling interaction which can be obtain by laser-induced- 
gauge-field method, A is the coupling coefficient, a are 



i.e., by the usual scalar product (-I-) combined with time- 
averaging. Hence, the quasienergies are obtained by 
computing the matrix elements of the Floquet operator 
H (t) — ihdt in the basis ([5]) with respect to the scalar 
product ([6]), and diagonalizing. By a straightforward 
calculation, we can obtain the matrix elements of some 
operators in Floquet Hamiltonian % (t): 



the Pauli matrices and d^ 

bond that connects site j to i. /i is the chemical po- 
tential and U < denotes an on-site attractive inter- 
action which is easy to obtain via an s-wave Feshbach 
resonance in cold atom system. The oscillating Hamil- 
tonian Hd, with K (t) = X (cos (ujt) , sin (ujt)), mimics a 
monochromatic electric dipole potential with frequency lo 
and amplitude K. This term can be realized experimen- 
tally by periodically shifting the position of a mirror em- 
ployed to generate the standing laser waves along x- and 
y-directions, and transforming to the comoving frame of 
reference [2J|. We choose i=l as the energy unit and 
the distance a between the nearest sites as the length 
unit throughout this paper. It was demonstrated in Ref. 
[3Qj that the Hamiltonian Hq in a mean field approxima- 
tion and combination with a perpendicular Zeeman field 
can support a topological superfluid. On the other hand, 
replaced the Hamiltonian Hq with an one-dimensional 
Bose-Hubbard Hamiltonian, it was shown in Ref. [23l 
that a phase transition from a superfluid to a Mott insu- 
lator can be induced by Hd in its one-dimensional form. 
When a Hamiltonian of quantum system has a peri- 
odic dependence on time, i.e., H {t)—H (t + T) with pe- 
riod T—2ti/uj, the Hamiltonian satisfies the discrete time 
translational symmetry, t ^ t + T, which can been de- 
scribed by Floquet's theorem (H, [sj, [Sij . Floquet's the- 
orem tell us that the Schrodinger equation with time- 
periodic dependent Hamiltonian has a complete set of 
solutions with the form !?/;„ (i)) = |u„ (t)) exp (— ie„i/?i). 
Here, the periodic function |u„ (t)) = |u„(t + T)), an 
analog of Bloch states known from spatially periodic crys- 
tals, satisfies the eigenvalue equation 



is a unit vector along the U{nj},m'\c\^Cja'\{ni\,m)\=e ''™ ™''^'' Jm'-m{zt]) ., 



[H{t)~ihdt]\un{t)) =e„|u„(t)) 



(4) 



We call Ti, {t)=H [t] — ihdt as the Floquet Hamiltonian 
and the eigenvalues £„ as quasienergies which are defined 
modulo the frequency aj=27r/T. 
The Floquet basis 



|{ni},m) = |{ni})exp 



hw 



di'K (t') • 2^ Vini+iraujt 
i 

(5) 



where |{nj}) indicates a Fock state with n, particles on 
the zth site, and m accounts for the zone structure [24| . 
consist of an extended Hilbert space of T-periodic func- 



{{n[},i 



\{ni} ,\ 



— TIIq^Q. 



m,m' 1 



(7) 
(8) 



({n-} , m'\ cJ^c|j^Q4,Ci-t. [{nj ,m)) ^ ni^niiSm^m', (9) 



where Jm'^m (zij) is 
(m' — m)th order and z, 



the Bessel function of the 



, K 



V 



Vt 
and taii9ij — 



IVi. 



^i] I yij 



(i"i)^ - (i"i)^, Uvj = {Yi)y - (rj) 

In the above matrix, the diagonal block of the Floquet 
Hamiltonian, 7^(™™)^ is the n-photon sector, i.e., the 
subspace with n photons and the non-diagonal blocks 
"Hi™ '"J with to' 7^ 771 correspond to the interaction be- 
tween different subspaces [27]. For sufficiently high fre- 
quencies, we can argue, from Eq. ((TUni), that the driven 
system ([T]) behaves similar to the undriven system ([5]), 
but with the tunneling matrix element t and the SO cou- 
pling A of the latter being replaced by the effective ma- 
trix element t ~ tJo (zij) and A ~ A Jo (zij), respectively. 
Now, suppose that we enhance the modulation strength 
K, then we have to consider the coupling of other pho- 
ton sectors. For simplify, we only consider coefficient 
of subspace with 77. = 1 photon on the subspace with 
77 = photon. When K is strong enough but still satisfy 



<< 1, the system has the effective Hamiltonian [2£ 



n 



eff 



^(00) 



1 

huj 



[n-',n+'] 



(10) 



Here, H^^ (^^^) denotes the non-diagonal block with 
mi—m — —1 (-1-1) around the 0-photon sector. According 
to Eqs. (IZlini), we can obtain 

^(00) = -f ^ Jo [z,,] c\c, -tXY,Jo (2.y) 4 (fT X d,,)^ r 



(y> 



(y> 



+ /i E 4^' + v^- E {4Ai + H.c.) , 

i i 

H-^ = -t5]e'^-J_i(z.,)c|c, 

iv) 

- iA^ e'*''^ J_i (zy) c| (a x d^ 



(11) 



ij <-■], 



(12) 






e-^^-J+i(z„)ct(axd,,) c,. (13) 



In the derivation of Eq. pT|) , we have made a mean field 
approximation and ips is the gap function. 

By using the Fourier transform of atomic operators c]^ , 
i.e., 






the Hamiltonian pOj) in square lattice system can be 
rewritten in the momentum space as 



^e//=^V',[(H.//(k)V'k, 



(15) 



where we have defined the four-component basis operator 
'0k=(ckt: Ckj,, c_j(.^, cLj^i )''". The effective Hamiltonian in 
(14) momentum space is given by 

I 



^e//(k) = 



£k-r(k,z) 2AJo(z)a(k) i^s \ 

2XJo{z)a*(k) ek + r(k,z) -t^^ 

-ips -e_k + r(k,.z) 2XJo{z)a*{k) 

^s 2AJo(z)a(k) -e_k-r(k,^)/ 



(16) 



r 



where a (k) 



sin fc„ 



i sin/cr, z 



method outlined in Ref. [30[, one can obtain a "dual" 



fs,r(k,z) = l^^J±2MI^li2l COS k, COS ky, and Hamiltonian 
£k = — 2tJo (z) (cosfc^ + cosfcj,) — /i. Following the 



H^ (k) = 



tA. - r (k, z) 2A Jo (z) a (k) 

2AJo(z)a*(k) -^, + r(k,z) Ek 

ek ^s+r (k, z) 

-£k -2A Jo (z) a (k) 



-£k \ 



-2A Jo (z) a* (k) 

-^, -r(k,z) / 



r 



(17) 



where the unitary transformation "H^ (k) — DH^ffD^ 
with 




It is easy to obtain the eigenvalues of Eq. (fT7|) with 

El — — v/i + f2^E2 
E3 
where we have defined 

A = r2 (k, z) + {el + i4) + 4J2 (z) A2 |a (k)!^ , 



I 2' -^2 — — V /l - /2, 

+V/i - /2; ^4 = +V/1 + /2, 



/2 = 2Vr2 (k, z) (4 + ^2) + 4AV2 (z) 4 |« (k)l^ 

It is obvious that if only tps ^ 0, i.e., the system lies 
in the superfluid phase, the energy levels Ei and E^, 
denoting the lowest and the highest band, will not touch 



each other. Next, we discuss the levels Ei and E2 (or E^ 
and E4). These two levels can touch each other if only 
the following two relations 



and 



r2(k,z) = o 



Jo (z) el (sin^ k^ + sin^ ky) = 0. 



are simultaneously satisfied. From F" (k, z) = 0, we have 
kx = ±§ or ky — ±-| in the first BZ. So, when Jo (z) 7^ 0, 
we have 

—fi = 2t Jo (z) (cos kx + cos ky ) . 

When —fi > |2iJo(z)|, the levels Ei and E2 
(or E3 and E4) will not touch each other for 
ever. Else, when — |2iJo(z)| < —fj, < |2iJo(z)|, 
El and £'2 (or E3 and E4) will touch each 

other at points ik^ ~ ±^,ky = arccos f 2tj^(z) ) ) ^^^ 
( krc = arccos ( 2tj^{z) ] ^^v ~ ^'Ir ^^ ^^^ following, we 



only consider the case, — /^ > |2iJo(2;)|, which means 
bands Ei and E^ will be off away from the other two 
levels E2 and E3, and then the topological properties of 
bands Ei and E^ will not change if we vary some pa- 
rameters. Therefore, we will only consider the topologi- 
cal properties of band E2 and E3, which have chances 
to contact each other at the high symmetry points, 
Ki=i,...,4 = (0,0);(0,7r);(7r,0);(7r,7r) in the first BZ, 
when satisfying the condition 



T^z)=s 



€ 



(18) 



because band-gap closing is an essential condition for 
the topological characteristic changes. We denote 
r2,^(K„z) = {l&\^Jl-y{z)/njjjf = V^iz). Consider- 
ing z << 1 and —ji > \2tJo {z)\, the two bands can only 



touch at point Ki = (0, 0) when varying the parameter 

z. 



III. TOPOLOGICAL PHASE TRANSITION 

We now study the topological properties of these two 
bands by two-band approximation at point Ki = (0, 0) in 
the first BZ. We will not consider the other points since 
the gaps at other high symmetry points will not shut 
down, leaving no influence to the topological changes. 
To have a basic idea of the topological features of the 
system, we explore it by using a two-band approximation 
at point Ki in the first BZ. We expand the Hamiltonian 
p?)) at point Ki and obtain 



n^, (q) = 



V's + r(z) 2XJo{z)q+ UJo{z) + ^i 

2AJo(z)g- -{^Js+T{z)) -{AtJo{z) + fi) 

-{4tJo{z)+n) iJs-r{z) -2AJo(z)g_ 

4tJo(z) + ^ -2XJo{z)q+ -{^Js-T{z)) 



(19) 



r 



with q±=qyzLiqx. Because f/'s < and F (z) > 0, when 
taking |4iJo {z) + fi\ <^ 1, we can see that the major 
contribution to bands £'2 and £'3 comes from the up- 
diagonal sector in above matrix and thus we may treat 
the others as a perturbation. Under this condition, we 
obtain an effective two-band Hamiltonian given by 

Hf/j.- (q) = 2AJo (z) qyax - 2AJo {z) q^a-y + M (z) a^, 

(20) 
where the corresponding mass term M {z)=^Js + T (z) + 
(4t Jo (z) + fi) I {ips — r (z)) and <yu=x,y,z the Pauli ma- 
trices. Let M {z)—0, we obtain the gapless condi- 
tion Eq. ([T8| again at Ki point. Eq. ([20|) can be 
written as 'HgA(q) — ad, where the vector d = 
{2XjQ{z)qy,—2XJo{z)qx,M{z)}. The topological fea- 
tures of the system can be characterized by the winding 
number ( first Chern number) of the Berry phase gauge 

field C = :^ / rffc^ / dfcj^d • (^ X ^) in the first Bril- 

louin zone, where d ~ d/|d|. When M {z) ^ 0, it is 
straightforward to obtain the winding number for the ef- 
fective system described by Eq. ((20)) . i.e., 



(21) 



C=-sign(Af(z)), 



This non-integral winding number appears since the devi- 
ations from this two-band approximation model at large 
momenta are not included in the above calculation of 
the winding number. So it can not be directly related 
to the topological features of the system; however, the 
change in the winding numbers is independent of the 
large- momentum contribution [36]. Let us discuss the 



change of the topological properties of superfluid system 
when we adjust the oscillating amplitude K of optical lat- 
tice. It is obvious that the initial non-driven system is in 
a trivial state which corresponds to z = and M (z) < 0. 
Now, we apply the driven field to the system and make 
M (z) < to M (z) > 0, the change in Chern number is: 

AC = i [sign {M (z)>o) - sign (M (z)^,)] = +1. (22) 

So, we get a topological superfiuid state with C — +1 for 
M (z) > . 

It is notable that gapless chiral edge states are usu- 
ally the hallmark of a topological system. Therefore, to 
further prove the above argument, we show the phase 
diagram and the band structures of the effective Hamil- 
tonian P^ in a striped geometry in Fig. 1. In Fig. 1(a), 
we plot the mass M (z) as a function of z, which can be 
adjusted by changing the modulation strength K. The 
point B denotes A-I (z) = with z = 0.1320, where the 
bands E2 and £3 contact each other at F point. The 
band structures in a striped geometry with 60 sites in x 
direction are shown in (b), (c) and (d) corresponding to 
the points A (z = 0), B (z = 0.132), and C (z = 0.2) 
in Fig. (a), respectively. In the numerical calculation, 
we take the typical parameters t = 1,A = 0.6, V's = 
—0.5, hoj — 0.05, and ^ = —4. It is clear that there is 
no edge state in the region / where the mass M (z) < 
with < z < 0.1320, so it is a trivial superfiuid. In 
contrast, there are a pair of edge states in the region // 
where the mass M (z) > with z > 0.1320, so it is a 
topological superfiuid. There is a phase transition from 
a trivial superfluid to a topological superfluid occurs at 





kyW 









0.1 0.2 

z=Ka/hco 


2 


/ 




HI 




-2 


'/"<\. 




FIG. 1: The phase diagram and the band structures of the 
system, (a) The mass M (z) as a function of z. The region / 
with (z) < is a trivial superfluid, while the region II with 
M (z) > is a topological superfluid. The band structures of 
the effective Hamiltonian pO|) in a striped geometry with 60 
sites in x direction are shown in (b), (c) and (d) corresponding 
to the points A (z = 0), B (z = 0.132), and C (z = 0.2) 
in (a), respectively. Other parameters f = 1, A = 0.6, i/'s = 
-0.5, &j = 0.05, and \i = -4. 



i/'s — ?► tpse^ , which is similar with the case in the topo- 
logical superconductor [38,]. Both cases give the similar 
Majorana fermion obviously confirmed from the Eq. (J16p 
and Eq. (|17p connected by the unitary transformation D. 
The zero mode solutions of the Majorana fermion can be 
obtained from the Bogoliubov-de Gennes (BdG) equa- 
tion, which has the similar form compared with that in 
Rcf. 30]. Moreover, such Majorana fermion excitations 
can be detected by the standard Raman spectroscopy 



IV. CONCLUSION 

In summary, we have discussed the topological super- 
fluid phase transition in a periodically driven square op- 
tical lattice. By using Floquet's theorem, we find that a 
Floquet topological superfluid will be created when the 
two-dimensional square optical lattice potentials are pe- 
riodically driven. This topological phase is interesting 
in hosting a Majorana fermion excitation which can be 
detected by Raman spectroscopy in cold atom system. 
Therefore we propose a novel scenario to create Majo- 
rana fermions which may play a key role in topological 
quantum computation. 



the point B where the gap is closed. 

Generally, there exists Majorana Fcrmionic excitation 
bounded with the vortex structure in the nontrivial topo- 
logical superfluid phase. Hence, we can obtain the 2-D 
Foquet Majorana fermions ^5| if there have the vortex 
structures in our system. The vortex structure can be 
produced from two different routes, one of which can 
be realized through the phase twist of the SO-produced 
lasers: A — > Ae*™^ with m the vorticity [30|- Another 
route is that the vortex structure can come from ini- 
tial rotation of the atomic cloud [33] ■ Then the vortex 
structure is coupled with the superfluid order parameter: 
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